Spinon excitation spectra of the J1J_1-J2J_2 chain from analytical calculations in the dimer basis and exact diagonalization

The excitation spectrum of the frustrated spin-$1/2$ Heisenberg chain is reexamined using variational and exact diagonalization calculations. We show that the overlap matrix of the short-range resonating valence bond states basis can be inverted which yields tractable equations for single and two spinons excitations. Older results are recovered and new ones, such as the bond-state dispersion relation and its size with momentum at the Majumdar-Ghosh point are found. In particular, this approach yields a gap opening at $J_2=0.25J_1$ and an onset of incommensurability in the dispersion relation at $J_2=9/17J_1$ [as in S. Brehmer \emph{et al.}, J. Phys.: Condens. Matter \textbf{10}, 1103 (1998)]. These analytical results provide a good support for the understanding of exact diagonalization spectra, assuming an independent spinons picture.Comment: 11 pages, 6 figure

Pair Density Waves in coupled doped two-leg Ladders

Motivated by Resonant X-ray scattering experiments in cuprate ladder materials showing charge order modulation of period $\lambda=3$ and 5 at specific hole densities, we investigate models involving the electronic t-J ladders and bosonic chains coupled via screened Coulomb repulsion. Extensive density matrix renormalization group calculations applied to the ladders/chains supplemented by a self-consistent mean-field treatment of the inter-ladder/chain coupling provide quantitative estimates of the charge order for $\lambda=3,4$ and 5. As previously proposed, such patterns correspond to the emergence of pair density waves which stem from the strong electronic correlations. We comment on the existence of a $\lambda=4$ modulation not seen so far in experiment.Comment: 4 pages, 4 figure

Magnetic responses of randomly depleted spin ladders

The magnetic responses of a spin-1/2 ladder doped with non-magnetic impurities are studied using various methods and including the regime where frustration induces incommensurability. Several improvements are made on the results of the seminal work of Sigrist and Furusaki [J. Phys. Soc. Jpn. 65, 2385 (1996)]. Deviations from the Brillouin magnetic curve due to interactions are also analyzed. First, the magnetic profile around a single impurity and effective interactions between impurities are analyzed within the bond-operator mean-field theory and compared to density-matrix renormalization group calculations. Then, the temperature behavior of the Curie constant is studied in details. At zero-temperature, we give doping-dependent corrections to the results of Sigrist and Furusaki on general bipartite lattice and compute exactly the distribution of ladder cluster due to chain breaking effects. Using exact diagonalization and quantum Monte-Carlo methods on the effective model, the temperature dependence of the Curie constant is compared to a random dimer model and a real-space renormalization group scenario. Next, the low-part of the magnetic curve corresponding to the contribution of impurities is computed using exact diagonalization. The random dimer model is shown to capture the bulk of the curve, accounting for the deviation from the Brillouin response. At zero-temperature, the effective model prediction agrees relatively well with density-matrix renormalization group calculations. Finite-temperature effects are displayed within the effective model and for large depleted ladder models using quantum Monte-Carlo simulations. In all, the effect of incommensurability does not display a strong qualitative effect on both the magnetic susceptibility and the magnetic curve. Consequences for experiments on the BiCu2PO6 compound and other spin-gapped materials are briefly discussed.Comment: 24 pages, 20 figure

Melting of a frustration-induced dimer crystal and incommensurability in the J_1-J_2 two-leg ladder

The phase diagram of an antiferromagnetic ladder with frustrating next-nearest neighbor couplings along the legs is determined using numerical methods (exact diagonalization and density-matrix renormalization group) supplemented by strong-coupling and mean-field analysis. Interestingly, this model displays remarkable features, bridging the physics of the J_1-J_2 chain and of the unfrustated ladder. The phase diagram as a function of the transverse coupling J_{\perp} and the frustration J_2 exhibits an Ising transition between a columnar phase of dimers and the usual rung-singlet phase of two-leg ladders. The transition is driven by resonating valence bond fluctuations in the singlet sector while the triplet spin gap remains finite across the transition. In addition, frustration brings incommensurability in the real-space spin correlation functions, the onset of which evolves smoothly from the J_1-J_2 chain value to zero in the large-J_{\perp} limit. The onset of incommensurability in the spin structure-factor and in the dispersion relation is also analyzed. The physics of the frustrated rung-singlet phase is well understood using perturbative expansions and mean-field theories in the large-J_{\perp} limit. Lastly, we discuss the effect of the non-trivial magnon dispersion relation on the thermodynamical properties of the system. The relation of this model and its physics to experimental observations on compounds which are currently investigated, such as BiCu_2PO_6, is eventually addressed.Comment: 13 pages, 13 figure

Slow quench dynamics of a trapped one-dimensional Bose gas confined to an optical lattice

We analyze the effect of a linear time-variation of the interaction strength on a trapped one-dimensional Bose gas confined to an optical lattice. The evolution of different observables such as the experimentally accessible onsite particle distribution are studied as a function of the ramp time using time-dependent exact diagonalization and density-matrix renormalization group techniques. We find that the dynamics of a trapped system typically display two regimes: for long ramp times, the dynamics are governed by density redistribution, while at short ramp times, local dynamics dominate as the evolution is identical to that of an homogeneous system. In the homogeneous limit, we also discuss the non-trivial scaling of the energy absorbed with the ramp time.Comment: 4 pages, 4 figures, version published in PR

Phase diagram of hard-core bosons on clean and disordered 2-leg ladders: Mott insulator - Luttinger liquid - Bose glass

One dimensional free-fermions and hard-core bosons are often considered to be equivalent. Indeed, when restricted to nearest-neighbor hopping on a chain the particles cannot exchange themselves, and therefore hardly experience their own statistics. Apart from the off-diagonal correlations which depends on the so-called Jordan-Wigner string, real-space observables are similar for free-fermions and hard-core bosons on a chain. Interestingly, by coupling only two chains, thus forming a two-leg ladder, particle exchange becomes allowed, and leads to a totally different physics between free-fermions and hard-core bosons. Using a combination of analytical (strong coupling, field theory, renormalization group) and numerical (quantum Monte Carlo, density-matrix renormalization group) approaches, we study the apparently simple but non-trivial model of hard-core bosons hopping in a two-leg ladder geometry. At half-filling, while a band insulator appears for fermions at large interchain hopping tperp >2t only, a Mott gap opens up for bosons as soon as tperp\neq0 through a Kosterlitz-Thouless transition. Away from half-filling, the situation is even more interesting since a gapless Luttinger liquid mode emerges in the symmetric sector with a non-trivial filling-dependent Luttinger parameter 1/2\leq Ks \leq 1. Consequences for experiments in cold atoms, spin ladders in a magnetic field, as well as disorder effects are discussed. In particular, a quantum phase transition is expected at finite disorder strength between a 1D superfluid and an insulating Bose glass phase.Comment: 24 pages, 23 figure

Statistical properties of the spectrum the extended Bose-Hubbard model

Motivated by the role that spectral properties play for the dynamical evolution of a quantum many-body system, we investigate the level spacing statistic of the extended Bose-Hubbard model. In particular, we focus on the distribution of the ratio of adjacent level spacings, useful at large interaction, to distinguish between chaotic and non-chaotic regimes. After revisiting the bare Bose-Hubbard model, we study the effect of two different perturbations: next-nearest neighbor hopping and nearest-neighbor interaction. The system size dependence is investigated together with the effect of the proximity to integrable points or lines. Lastly, we discuss the consequences of a cutoff in the number of onsite bosons onto the level statistics.Comment: 18 pages, 15 figure

Reply to "Comment on `Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined' ''

In his Comment [see preceding Comment, Phys. Rev. A 82, 037601 (2010)] on the paper by Roux [Phys. Rev. A 79, 021608(R) (2009)], Rigol argued that the energy distribution after a quench is not related to standard statistical ensembles and cannot explain thermalization. The latter is proposed to stem from what he calls the eigenstate thermalization hypothesis and which boils down to the fact that simple observables are expected to be smooth functions of the energy. In this Reply, we show that there is no contradiction or confusion between the observations and discussions of Roux and the expected thermalization scenario discussed by Rigol. In addition, we emphasize a few other important aspects, in particular the definition of temperature and the equivalence of ensemble, which are much more difficult to show numerically even though we believe they are essential to the discussion of thermalization. These remarks could be of interest to people interested in the interpretation of the data obtained on finite-size systems.Comment: 3 page